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We show that dissipative quantum state preparation processes can be protected against qubit 
dephasing by interlacing the state preparation control with dynamical decoupling (DD) control 
consisting of a sequence of short vr-pulses. The inhomogeneous broadening can be suppressed to 
second order of the pulse interval, and the protection efficiency is nearly independent of the pulse 
sequence but determined by the average interval between pulses. The DD protection is numerically 
tested and found to be efficient against inhomogeneous dephasing on two exemplary dissipative 
state preparation schemes that use collective pumping to realize many-body singlets and linear 
cluster states respectively. Numerical simulation also shows that the state preparation can be 
efficiently protected by vr-pulses with completely random arrival time. Our results make possible 
the application of these state preparation schemes in inhomogeneously broadened systems. DD 
protection of state preparation against dynamical noises is also discussed using the example of 
Gaussian noise with a semiclasscial description. 
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I. INTRODUCTION 

Dissipative quantum state preparation has recently 
emerged as a conceptually new approach for realizing re- 
sources of multipartite quantum entanglement. In the 
control of an open quantum system, dissipative chan- 
nels are usually the sources of errors and play delete- 
rious roles. However, by properly tailoring the dissipa- 
tive channels in an open quantum system, the irreversible 
dynamics can have steady states containing the impor- 
tant resource of quantum entanglement. Earlier studies 
in small scale systems have shown that two-qubit entan- 
glement can be generated in the steady states of various 
dissipative dynamics driven either by external incoherent 
pumping or by the tailored coupling with the environ- 
ment [iH^. Recent studies have lead to the discovery of 
various dissipative state preparation approaches for real- 
izing multipartite entanglement in large scale systems (Bi- 
ll 3]. In particular, it has been shown that by using only a 
few incoherent pumping channels that couple collectively 
to the qubits, various types of multipartite entanglement 
can be prepared, including the linear cluster states [lO|, 
many-body singlet states [11], and symmetric and asym- 
metric Dicke states 0,ll3]. These multipartite entangled 
states can be the crucial resource for measurement based 
quantum information processing fl3-[T7|. With the de- 
sired quantum states unconditionally realizable as steady 
states of the irreversible dynamics, the dissipative state 
preparation schemes have the advantage of robustness 
and do not need accurate temporal controls as compared 
to conventional state preparation by coherent evolutions. 

Most dissipative state preparation schemes are pro- 

Sosed for qubit systems with homogeneous resonances [1- 
, [iqI-[i3|. In reality, the qubit resonance can have in- 
evitable inhomogeneous broadening which leads to in- 
homogeneous dephasing. Moreover, fluctuation of the 
environment can also induce dynamical noise that leads 
to dephasing as well. These unwanted noise channels will 
compete with the tailored dissipative channels for prepar- 



ing quantum entanglement in the irreversible dynamics. 
The fidelity of the steady state with the target quantum 
state will decrease with the increase of the noise strength 
and the number of qubits. These inevitable noises will 
render the state preparation scheme impractical in large 
scale systems in the presence of inhomogeneous broad- 
ening and dynamical noise. Whether or not we can sup- 
press the noises while preserving the desired dissipative 
channels in the irreversible dynamics becomes the key. 

Dynamical decoupling (DD) has been an extensively 
explored approach in high-precision magnetic resonance 
spectroscopy [18-21], and for suppressing the effect of 
qubit dephasing in quantum information processing [22\- 
25]. The basic idea is to introduce a sequence of pulses 
on the qubit in order to average out the inhomogeneous 
broadening and the unwanted coupling with the environ- 
ment, hence eliminating the noise effects on the qubit. 
Past studies on dynamical decoupling have focused on 
how to freeze the evolution of a qubit prepared on an 
arbitrary quantum state such that the quantum memory 
time can be enhanced. One efficient DD approach for 
protecting quantum memory uses sequences of instanta- 
neous TT-pulses with properly chosen arrival time. Inho- 
mogeneous dephasing can be removed transiently at a 
certain time known as the spin echo time, while dynam- 
ical noises can also be suppressed to high orders at the 
spin echo times depending on the pulse sequences. The 
protection of coherent evolutions by DD control pulses 
has also been studied l26l-l3Ql. 



In this paper, we show that dissipative quantum state 
preparation processes can be protected against qubit de- 
phasing by interlacing the dissipative control with DD 
control consisting of sequences of 7r-pulses. The DD 
control is introduced here to suppress the qubit de- 
phasing while preserving a nontrivial irreversible evo- 
lution for generating entanglement in the steady state. 
For suppression of inhomogeneous dephasing, we com- 
pare DD controls consisting of repetitions of a basic 
unit which takes various designs including the Carr- 
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Purcell-Meiboom-Gill (CPMG) (HI, HI, concatenated 
(CDD) [33-^8], and Uhrig sequences (UDD) ^^^l We 
found the leading noise term in the Magnus expansion 
is of second order of the pulse interval and the residue 
effect of inhomogeneous dephasing is largely determined 
by the average pulse interval only, nearly independent 
of the types of the DD pulse sequences. High fidelity 
state preparation is therefore possible even in the pres- 
ence of substantial inhomogeneous broadening when the 
DD pulse sequence with sufficiently small pulse interval is 
applied. DD protection from dynamical noise is discussed 
for the example of Gaussian noise of a semiclassical de- 
scription, where the leading noise effect is also found to 
be of second order of the pulse interval. The order of 
noise suppression is qualitatively different from that in 
the DD protection of quantum memory because of the 
presence of the state-preparation channel, and its inter- 
ference with the noise channels determines the residue 
noise effect. 

For the scheme of preparing the many-body sin- 
glet [nl, [lH, we perform systematic numerical simula- 
tions of the DD protection against inhomogeneous de- 
phasing using CPMG, concatenated, and Uhrig pulse se- 
quences. We find that the Magnus expansion converges 
well and the residue inhomogeneous broadening effect is 
well accounted by the leading noise term in the Magnus 
series. Numerical simulation confirms that the protection 
efficiency is nearly independent of the pulse sequence and 
depends only on the average pulse interval. Moreover, the 
application of 7r-pulses with completely random arrival 
time (43I , [43 | is found to be efficient in suppressing the 
effect of inhomogeneous broadening. Thus, the DD pro- 
tection is robust against errors in the pulse arrival time. 
We also explain how to implement the DD protection on 
the scheme of preparing linear cluster states [10], where 
our numerical simulation confirms that high fidelity state 
preparation can be realized in the presence of substan- 
tial inhomogeneous broadening through the DD control. 
These dissipative preparation schemes for preparing the 
two important classes of multipartite entanglement are 
therefore applicable to inhomogeneously broadened qubit 
systems. 

The paper is organized as follows. In section II, we 
discuss the basic idea of protecting dissipative quantum 
state preparation against noise channels by interlacing 
the state-preparation control with the DD control pulses. 
The evolution of the system is expanded in terms of the 
Magnus series. In section III, we discuss the DD pro- 
tection of the state preparation against inhomogeneous 
dephasing, and derive coefficients of Magnus series for 
investigating the order of noise suppression. In section 

IV, we discuss the DD protection on the preparation of 
many-body singlet states against inhomogeneous dephas- 
ing and show numerical simulation results. In section 

V, we numerically study the DD protection of the lin- 
ear cluster states in the presence of inhomogeneous de- 
phasing. In section VI, DD protection against dynamical 
noise is discussed for the example of Gaussian noise of 



a semiclassical description. Finally, the conclusions are 
summarized in section VII. 



II. MAGNUS EXPANSION 

In the dissipative quantum state preparation, the evo- 
lution of the open quantum system is described by the 
master equation as p = >Cp [p], where p is the density 
matrix of the system and Cp [p] is the super-operator 
describing the dissipative dynamics tailored for generat- 
ing entanglement. In a t^ical dissipative quantum state 
preparation process |ll-[l3|, the dissipative channels Cp 
are time independent and are engineered such that the 
desired entangled states are obtained as the steady states 
of the irreversible dynamics. 

In reality, there exist other noise channels which also 
affect the system dynamics as denoted by the super- 
operator £n [p], which can have various origins such as 
inhomogeneous broadening of the level spacing, and dy- 
namical ffuctuation from the coupling with an environ- 
ment. The unwanted evolution £n [p] will compete with 
the tailored dissipative dynamics Cp [p]. As a result, the 
steady state of the dynamics p = Cp [p] + £n [p] deviates 
from the target entangled state. The state preparation 
fidelity will drop with the increase of the noise strength 
and the number of qubit s. To achieve high fidelity state 
preparation in large scale systems, the undesired noise 
channels C^ need to be suppressed while preserving the 
desired dissipative dynamics Cp [p] at the same time. 

We consider interlacing the state-preparation control 
with the DD control consisting of a sequence of short 
pulses which realize the unitary rotations Ui at time 
on all qubits. If the duration of the pulse is much smaller 
than the pulse interval, the rotation can be considered 
as instantaneous. It is convenient to analyze the evo- 
lution of the system in the toggling frame that follows 
these rotations of the qubits. In the toggling frame, the 
super-operators for the noise channels and the desired 
dissipative channels become: C^ [t; p] = U{t)C-^ [p] U{ty 
and £T [^.p] = U{t)Cp [p] /7(t)t, where U{t) = Uj 
for Ti < t < Ti+i. Below we focus on DD control using 
sequences of short tt pulses (see Fig. [1]). In such case, 
Ui = (Jx^ and we have U{t) = Gx in the odd intervals 
(i.e. Ti <t < Ti+i with odd z), and U{t) = I in the even 
intervals (i.e. Ti < t < r^+i with even i). 

If the noise super-operators in the toggling frame have 
overall sign changes for adjacent intervals as C^ [t] p] = 
f{t)C^[p], where f{t) = {-!)' for n < t < r^+i, the 
effect of the noise can then be averaged out over time. 
This is the case for pure dephasing noise. However, in 
the protection of dissipative quantum state preparation, 
the desired dissipative channels for generating entangle- 
ment shall also be preserved. This imposes an additional 
requirement [26j, i28j. The simplest scenario is to require 
the dissipative channel to always take the desired form 
Cp^o in the toggling frame, i.e. Cp [t; p] = Cp^o [p]. Thus, 
in the laboratory frame, the super-operator for the dis- 
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FIG. 1. (Color online) Schematics of interlacing DD con- 
trols with the dissipative quantum state preparations. In the 
toggling frame, the noise has a sign change between odd and 
even time intervals and the noise effects average out over time. 
The effects of the dissipative preparation channel in the odd 
and even time intervals constructively add. In the laboratory 
frame, the dissipative preparation channel shall then take two 
different forms for the odd and even time intervals in general. 



sipative channel shall alternate between the two forms 
>Cp,o [p] and £p q [p] = (TxJC.p^o [p] cf^ in the even and odd 
intervals respectively (see Fig. [1]). In general, i2p^o and 
C'y> are different operators. Thus, for a dissipative state 
preparation protocol to be protectable by DD controls 
consisting of tt pulses, the dissipative dynamics shall 
switch between the two forms each time a tt pulse is ap- 
plied if [cFx^ £p [p]] ^ 0. This can indeed be implemented 
in most dissipative preparation protocols as the dissipa- 
tive channels Cy> [p] are engineered ones. 

In the toggling frame, the super-operator for the dis- 
sipative preparation channel is time-independent, hence 
the effect adds constructively over time. The super- 
operator for the noise channels has alternating signs for 
even and odd time intervals, hence the effect destruc- 
tively cancels. The master equation in the toggling frame 
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FIG. 2. (Color online) Schematics of DD control, formed by 
repetitions of a basic unit with duration tp. The basic unit 
consists of a A/'-pulse sequence which can take various designs, 
including the CPMG, concatenated and Uhrig sequences. 



three terms of this series are: 



Jo 



t ptl pt 



03^ = ^1 I I {[C[h],[c[t2],CM 

[C [t3],[C[t2],C[h]]]}dhdt2dt3, 



(4a) 
(4b) 

(4c) 



where C[t] = >Cp,o + /(^)>Cn, and [C[ti] , C[t2]] p = 
jC [ti; C [t2; p]] — jC [t2; jC [ti; p]] is the Poisson bracket for 
the super-operator £ [t] . The Magnus expansion is con- 
vergent if /q \\jC \t'] II dt' < TT, where ||£|| is a two- form of 
super-operator C [45|. It is convenient to investigate the 
suppression of noise channels with the Magnus expan- 
sion. The residue effects of noises are characterized by 
the leading term in the Magnus series that contains the 
noise operators, as we will discuss explicitly in the next 
section for inhomogeneous dephasing noise. 



/5 = £p,o[p] + /(t)£N [p], 



(1) 



which has the formal solution: 



III. PROTECTION OF DISSIPATIVE STATE 
PREPARATION FROM INHOMOGENEOUS 
DEPHASING 



p (t) — Texp 



l\cp^o^f{t')C^)dt' 
Jo 



p{0) (2) 



with time ordering T and initial density matrix p (0). Eq. 
(IH) can be rewritten in terms of the Magnus expansion: 



(3) 



Here, the formal evolution operator is expanded in terms 
of the Magnus series as Q {t) = XlfcLi (^), and the first 



In this section, we investigate the protection of dissi- 
pative state preparation from inhomogeneous dephasing 
by the DD controls. The super-operator for this noise 
channel is 



>Cn [p] 



ujicr^ , p 



(5) 



where the inhomogeneous dephasing is the consequence 
of inhomogeneous broadening in the qubit resonance uOi. 
The first three terms of the Magnus expansion then be- 
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come: 

(t) = t [Cp,o + ci (t) £n] , (6a) 
!^2W = t'c2(t)[£p,o,£N], (6b) 

1^3 (t) = t^C3,a (t) [Cp,o, [Cp ,0, >Cn]] 

+t^c3,6(t)[/:N,[/:p,o,/:N]], (6c) 

where [>Cp,o, >^n] P = jCp,o [>^n [p]] - >^n [>^p,o [p]] stands 
for the Possion bracket for the super-operators. Because 
of the interference between the super-operators £p,o and 
£n, these Poisson brackets are typicahy nonzero. 



The coefficients in Eq. (|6]) are given by: 

1 r' 

ci{t) = - / f{ti)dh, (7a) 
^ Jo 

C2(t) = ^ / ^^3{ti)dti, (7b) 
1 

C3,a(t) = tlh{h)dh, (7C) 

C3,6 W = / ^? [6^2 (tl) / (ti) - Ci (ti) ^ (ti)] , 

(7d) 

Here^(ti) = ci (ti)-/(ti) and /i(ti) = 6c2 (ti)-^(ti). 
Since /(ti) = (—1)* for ti < ti < r^+i, if the total du- 
ration of the odd intervals equals that of the even inter- 
vals, ci will vanish and the leading term of the Magnus 
expansion Vti only contains the ideal dissipative prepa- 
ration channel we wish to preserve. The residue effect 
of the inhomogeneous broadening is given by the higher 
order Magnus terms, which results from the interference 
of the noise channel with the preparation channel. 

In the dissipative state preparation, the desired quan- 
tum state is obtained as the steady state of the dynam- 
ics, and the timescale to reach the steady state is de- 
termined by the dissipative channel i2p,o- We consider 
a practical scenario of interlacing the DD control with 
the dissipative preparation (see Fig. [2]). The DD control 
has a periodic structure which consists of repetitions of a 
basic unit with the duration tp. The basic unit is an N- 
pulse sequence that can take various designs, including 
the CPMG (3i|,[32|, concatenated ^33.-^1, and Uhrig [39- 
|42| sequences. The basic unit of the DD control can be re- 
peated as many times as necessary along with the dissipa- 
tive evolution until a steady state is reached. The average 
interval between the 7r-pulses can be defined f = tp/N. 

For a periodic control with period tp, it has been 
proved that Q^k (Itp) = l^k {tp)^ where / = 1,2,.... (45| . 
Therefore the coefficients Cn possess the scaling law 

c„ {Hp) = ;i-"c„ {tp) = (^-f^ an. (8) 

Here an = Cn{tp) depends on the structure of the A/'-pulse 
unit, but is independent of its duration tp as is evident 



from Eq. ([7j). The Magnus terms then have the form 

^1 {t = Hp) = t [jCp^o + aiC^] , (9a) 

^2 {t = Up) = t{a2tp) [/:p,o, jC^] , (9b) 
Qs {t = Hp) = t {as^atl) [>^p,o, [>^p,o, >Cn]] 

+^ {(^sMI) [^n, [>Cp,o, jCn]] . (9c) 

Obviously, for the system evolution in a given duration 
the unwanted dynamics described by the higher order 
Magnus terms ^k>2 can be better suppressed when tp is 
smaller for a given type of A/'-pulse unit. This is equiv- 
alent to using pulse sequences with shorter interval f. 
When tp is sufficiently small, the residue effects of noises 
are predominantly captured by the leading noise term, 
i.e. ftk>2 with the smallest /c, and convergence of the 
Magnus expansion can be expected. 

In table I, we list the coefficients an for the first sev- 
eral Magnus terms, where the basic unit uses various DD 
pulse sequences. In the DD protection of quantum mem- 
ory, it is well known that inhomogeneous dephasing can 
be completely removed, and dynamic noises can be sup- 
pressed to arbitrarily specified order of the pulse interval 
by using higher order concatenation design [33-38J or the 
Uhrig design [39-42]. Here for the protection of the state 
preparation, we find a qualitative different behavior. Be- 
cause of the presence of nontrivial dissipative dynamics 
for the state preparation, the inhomogeneous broadening 
always has residue effects due to its interference with the 
preparation channel (c.f. Eq. ([6])). For all pulse sequences 
considered, ^3^5 can only be suppressed but never vanish, 
and the inhomogeneous noise can be suppressed at most 
to the second order of the pulse interval. 

An interesting observation is as^^A^^ is nearly inde- 
pendent of the type of pulse sequences. We note that the 
Magnus expansion in such cases reads 

n{t) = t {Cp^o + as.biV'r^ [£n, [£p,o, -Cn]] + 0(f 3)) 

(10) 

where we have used tp = Nr. Thus, the efficiency on 
the suppression of inhomogeneous dephasing is mostly 
determined by the average pulse density f, instead of the 
order of the CDD or UDD sequences. 

Here we note that Cp^q in the above derivation can 
stands for a general evolution that one wishes to preserve 
while using DD control to suppress the dephasing. 

IV. DD PROTECTION ON THE DISSIPATIVE 
PREPARATION OF MANY-BODY SINGLETS 

In this section, we numerically study the DD pro- 
tection on an exemplary dissipative state preparation 
scheme in the presence of inhomogeneous dephasing. The 
scheme uses collective pumping to prepare many-body 
singlets flT]. It has been shown that when an ensem- 
ble of spins are collectively pumped by the homogeneous 
collective lowering operator s~ and collective raising 
operators ^n^^i^n with inhomogeneous coefficients a^, 
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TABLE I. Coefficients of leading Magnus terms in Eq. Q. The coefficients are evaluated for one basic unit of evolution with 
duration tp where the DD control is a A^-pulse sequence taking various designs (see text) . 
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FIG. 3. (Color online) Steady state population of many-body 
singlet P( J = 0) versus the strength of inhomogeneous broad- 
ening A, when the preparation process is protected using dif- 
ferent DD pulse sequences. The ensemble consists of 6 qubits, 
with inhomogeneously broadened energies Ui = A(7 — 2i)/5, 
i = 1,2, ...,6. The collective pumping rates are chosen as 
= 10A^ (see text). The red solid, blue dashed, brown 
dotted, orange dashed, and light yellow dotted curves show 
numerical results by solving the exact master equation, where 
the basic unit of the DD control consists of the CPMG, CDD3, 
CDD4, UDD5 and UDDIO pulse sequences respectively. For 
the simulation presented in (a), the duration of the basic unit 
tp = 10~^A~^ for all DD controls compared, as shown in (c). 
In (b), the duration of the basic unit tp — 10~^A~^A^, N be- 
ing the number of pulses in the basic unit, and the average 
pulse interval r = tp/N is taken to be the same for all DD 
controls compared, as shown in (d). 



the steady state is a highly entangled one in the neigh- 
borhood of many-body singlets. In the absence of inho- 
mogeneous and homogeneous dephasing, the singlet pop- 
ulation in the steady state P(J = 0) > 20%. Here we 
consider the scenario on the preparation of many-body 
singlet as investigated in Ref. [12]. An spin ensemble 
with the even number of spins is initially in the fully po- 
larized state. The only two collective pumping operators 
needed for preparing the many-body singlet from this 
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FIG. 4. (Color online) Steady state population of many-body 
singlet P( J = 0) as function of the inhomogeneous broaden- 
ing A and the average pulse density ft — N/tp, where the 
basic unit of DD control uses various A/'-pulse sequences. The 
ensemble consists of 6 qubits, with inhomogeneously broad- 
ened energies uji = A(7 — 2i)/5, i = 1, 2, . . . , 6. The collective 
pumping rates are chosen as A^ = lOOA^. 



initial state are respectively s~ and Xln(~)^'^n • "^^^ 
preparation channel Cp^o is simply realized by the simul- 
taneous pumping with these two operators, described by 
the Lindblad form: 



j=l,2 



L,pL] 



(11) 



The two Lindblad operators are Li = y/A^^^^n 
1/2 = V^^n{—)^s^^ where Ah and A^ characterize the 
strength of the pumping. In the absence of dephasing 
noises, the population of many-body singlets P{J = 0) ^ 
40% in the steady state when Ah ^ is satisfied. Here, 
A~^ characterizes the timescale for reaching the steady 
state. 

In the presence of inhomogeneous broadening in the 
qubit resonance, the steady state will deviate from the 
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FIG. 5. (Color online) Convergence of the Magnus expansic 
Solid and dot-dashed curves show the steady state populati 
of many-body singlet P(J = 0) obtained by exactly solvi 
the master equation. Dashed and dotted curves are stea 
state P(J = 0) obtained from the Magnus expansion Eq. Q 
where we only keep the leading higher order Magnus tei 
(see text). The ensemble consists of 6 qubits, with inhomo^ 
neously broadened energies uoi = A(7 — 2i)/5, i = 1, 2, ... , 
The collective pumping rates are chosen as Ah = lOOAi. T 
blue, orange, red and green solid lines represent numeric 
results when the basic unit of DD control uses the CPM 
UDD3, CDD6 and UDD42 pulse sequences, respectively. T 
average pulse density is chosen as n = N/tp = lO^'^^A^. 



target one because of the inhomogeneous dephasing. \ 
use the population of singlets P(J = 0) in the stea( 
state as the figure of merit. To preserve the collect! 
pumping channel while applying the DD control, in t 
laboratory frame, the Lindblad operators shall alterna 
between two forms in the even and odd intervals of t 
DD pulse sequence: 



^Y.n^n^ even interval, 
Zln interval. 



and 



L2 (t) = 



/AlEn(-)"s; 



even interval, 
odd interval. 



(12) 



(13) 



The numerical simulations are shown in Fig. 3-6. 
In Fig. [21 we show the state preparation efficiency with 
the increase of the inhomogeneous broadening when the 
preparation process is protected with various pulse se- 
quences. The steady state population of many-body sin- 
glet P( J = 0) is obtained by solving the master equation 
exactly. From Fig. [3](a), apparently, the state prepara- 
tion can tolerate larger inhomogeneous broadening when 
the basic unit of DD control uses higher order concatena- 
tion designs or Uhrig designs when tp is fixed. Interest- 
ingly, nearly identical protection efficiencies are observed 
for CDD3 and UDD5, and also for CDD4 and UDDIO. 
We find the common feature between CDD3 and UDD5, 
and between CDD4 and UDDIO is that they consist of 



FIG. 6. (Color onUne) Steady state population of the many- 
body singlet P(J = 0) as function of the inhomogeneous 
broadening A and the average pulse density n — N/tp when 
the DD control uses completely random pulse sequence. The 
ensemble consists of 6 qubits, with inhomogeneously broad- 
ened energies uji = A(7 — 2i)/5, i = 1, 2, . . . , 6. The collective 
pumping rates are chosen as = lOOAi. 



the same number of tt pulses. This is consistent with the 
analysis in section HIH namely, the efficiency of protection 
from inhomogeneous dephasing is largely determined by 
the density of pulses, rather than the actual design of the 
sequence. In Fig. [3]^c), we compare the performance of 
various DD controls by fixing the average pulse interval 
f = tp/N^ and different pulse sequences indeed exhibit 
nearly identical protection efficiency. 

In Fig. m we examined the protection efficiency as 
function of both the inhomogeneous broadening and the 
average pulse density n = 1/f. For the four DD con- 
trols compared, nearly identical protection efficiency is 
seen. Moreover, in the double log^^g plot, the contour 
lines are nearly straight lines with slope 1, showing that 
the strength of the residue noise scales as A^f^. This sug- 
gests that as^i)N'^f'^ [>Cn, [>Cp,o, ^n]] is indeed dominating 
over other higher order Magnus terms and captures the 
major effect of the residue noise (c.f. Eq. (p!Q|) ). 

In Fig. [SI we check the convergence of the Magnus ex- 
pansion by comparing the steady state P( J = 0) numer- 
ically solved from the exact master equation Eq. ([T]) with 
that obtained from the Magnus expansion Eq. (|3]) where 
we only keep the leading higher order Magnus term, i.e. 
by taking n{t) = t {Cp^o ^ as^bN^f^C^.lCp^oXN]]) • 
Excellent convergence is found even when the inhomo- 
geneous noise is strong enough to diminish the steady 
state population of the singlets. 

Motivated by the finding that the protection efficiency 
is determined by the average pulse density and is nearly 
independent of the pulse structure, we further tested the 
possibility of protecting the state preparation process by 
a completely random sequence of 7r-pulses. Since the DD 
control in such case does not have the periodic structure, 
the analytical results of Eq. ^ and Eq. ([9|) are no longer 
applicable. Nevertheless, as shown in Fig. [6l the state 
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preparation can still be well protected from inhomoge- 
neous dephasing. The slope of the contour line is close to 
2 in the shown parameter regime, from which we can tell 
that the residue noise scales as A^f, different from the 
periodic DD controls where the basic unit uses CPMG, 
CDD and UDD designs. 

Our numerical study shows that sequences of tt pulses 
with random arrival times are also efficient in protecting 
the dissipative state preparation. This means that the 
protection is insensitive to the errors in the arrival time 
of the TT pulses and accurate control on the arrival times is 
not necessary. The only temporal control needed is that 
the TT pulses have to be synchronized with the switches 
of the pumping operators in the laboratory frame. 



V. DD PROTECTION ON THE DISSIPATIVE 
PREPARATION OF LINEAR CLUSTER STATES 

It has also been shown that collective pumping can be 
used to prepare linear cluster states in an atomic ensem- 
ble embedded in a cavity tlO|. Here we investigate the 
DD protection of this dissipative preparation scheme. 

Key to the preparation of linear cluster states is to re- 
alize the stabilizers through the competition between the 
optical pumping and the atomic spontaneous emissions. 
The level scheme for realizing the stabilizer Ss = ^^f spg 
is shown schematically in Fig. [71 The qubit is defined 
in the ground state manifold {\L)j^ \R)j} of a three level 
atom. Optical fields denoted by the double head solid 
lines in Fig. [71 couple the two ground states to the com- 
mon excited state \E)j with Rabi frequencies and 
respectively, where j labels the atoms. All atoms are cou- 
pled to a common cavity mode with coupling strength 
g. and are the spontaneous emission rate from 
\E)j to the two ground states respectively, whose values 
are usually considered as the same. Those spontaneous 
emissions are denoted by wavy lines in Fig. [71 When the 
driving fields are tuned in resonances with the transitions 
as indicated in the figure, the eigenstate of stabilizer Ss 
is realized as the steady state of the pumping. Together 
with single qubit rotations, a complete group of stabiliz- 
ers for n— qubit linear cluster states can be realized one 
by one (e.g., Sl = sfs^, Sj = sfsfsf, = spfsl and 
Sf = S3 54 for 4-qubit linear cluster states). And thus 
the fidelity of the linear cluster states can approach unity 
when damping rate of atomic spontaneous emission F is 
much smaller than the cavity- atom coupling strength g. 
More details of the preparation scheme can be found in 
Ref. fw]. 

The scheme relies on delicate engineering of detun- 
ings, hence the preparation fidelity will be substantially 
affected by inhomogeneous broadening of the qubit reso- 
nances. We consider interlacing DD pulse sequences with 
the pumping to suppress the effect of this noise. Since 
the TT— pulse on the qubit simply switches \L)j and \R)j 
as to preserve the desired form of pumping channel in 
the toggling frame, one simply needs to switch values 
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FIG. 7. (Color online) Level scheme for realizing the sta- 
bilizer S3 = In the notations such as \RLLO), the 
first three letters denote the state of qubit 1, 2, 3 respec- 
tively and the last integer denotes the number of cavity pho- 
tons. The four involved states in the excited state manifold 
are: |1) = 1/V2{\RLL1) + \ELL0)), |2) = 1/V2\RLR1) - 
1/2{\ELR0) + \RLEO)), |3) = 1/V2\RRL1) - 1/2{\ERL0) + 
\RELO)), |4) = 1/V2\RRR1) - 1/V6{\ERR0) + \RERO) + 
\RREO)). The wavy lines denote the spontaneous emission. 
Optical fields resonantly pump the four transitions denoted 
by the double head solid lines. The eigenstate of stabilizer 
5*3 is realized as the steady state of the pumping, (a) and 
(b) correspond respectively to pumping controls in the labo- 
ratory frame before and after a tt pulse is applied. In (a) the 
cavity couples to the \R) ^ \E) transitions with strength g^, 
while in (b) the cavity couples to the the \L) ^ \E) transi- 
tions with strength g^ [46| . In the toggling frame, the two 
pumping controls realize the same dissipative channel. 



for every pair of Rabi frequencies ftj" and ftj^ and corre- 
sponding frequencies of the pump fields in the laboratory 
frame whenever a 7r-pulse is applied in the DD control. 
The required switching of the controls in the laboratory 
frame is illustrated in Fig. [71 where part (a) and (b) cor- 
respond respectively to the control before and after a 
TT-pulse is applied. This is equivalent to using different 
optical pumping fields in the even and odd intervals of 
the DD control. 

In Fig. [HI we show numerical simulation on the pro- 
tection of preparing a 4-qubit linear cluster state. The 
parameters = = O.Olg and nf = nf = O.Olg. 
The basic unit of the periodic DD control uses the UDD3 
sequence. Indeed, the DD protection is efficient in pro- 
tecting this dissipative preparation scheme when the in- 
homogeneous broadening A < 0.1^. The fidelity of the 
steady state with the target one approaches unity when 
the average pulse interval f is small. 
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Fidelity 




FIG. 8. (Color online) Fidelity between the steady state and 
the target 4-qubit linear cluster state under the DD protection 
where the basic unit uses the UDD3 sequence. The inhomo- 
geneously broadened qubit resonances are Ui — A(5 — 2i)/3, 
I — 1, 2, 3, 4. The parameters for the preparation channel are: 
VL^ = Vtf = 0.01^, = rf = 0.01^. See text. 



the DD control used jH, 

\F {UJ, t)\ = E / P ^^fitl)fitl-t2) COSUt2C {t2) . 

(17) 

In the case of protecting quantum memory, the prepa- 
ration channel >Cp,o is absent. Thus O (^2) = (>Cy) , 
which is time-independent and can be taken out of the 
integral in Eq. ([T7j). In such case, concatenated DD or 
Uhrig DD pulse sequences can generate the filter function 
that suppresses dynamic noise to arbitrarily high order 
of pulse intervals [39, 48- 50]. 

However, for the protection of a dissipative state prepa- 
ration, because of the presence of the nontrivial prepara- 
tion channel >Cp,o, (^2) can be time-dependent due to 
the interference between Cp^o and Cy. The effect of dy- 
namics noises can not be arbitrarily suppressed as in the 
case of DD protection of quantum memory. In general, F 
scales as f^, and the dynamical noise can be suppressed 
to the second-order of the pulse interval. 



VI. DD PROTECTION FROM DYNAMIC 
NOISES 

Dynamic fluctuation in the environment can also lead 
to qubit dephasing. Here we investigate the possibil- 
ity of DD protection of the dissipative state prepara- 
tion in the presence of dynamic noise, for which we 
assume a semiclassical form. The qubit energies are 
given by Hsc = ^iBi{t) (jf ^ where Bi{t) is stochas- 
tic and is characterized by the autocorrelation functions 
{Bi (t) Bj {t')) = G{t- t') 6ij. We further assume the 
noises are Gaussian, for which we can derive the Kubo 
stochastic Liouville equation [47] in the toggling frame as 

p=Cp,o[p] + Cf[t;p] (14) 

where the effective noise operator is 

. Jo 

(15) 

with O (t) = CyCy{t), where Cy[p] = —i[af,p] and 

£VW = exp(/:p,ot)/:Vexp(-/:p,ot). fit) = i-y for 
Ti < t < Ti+i is the step-like function corresponding to 
the DD control. 

Following Eq. (|4a|) . the first term of the Magnus series 
can be written as 

1^1 (t) = t£p,o [p] + / G{u)\F{u,t)\dcu, (16) 

where G {uj) = l/27r J^^G (t) exjp{—iujt)dt is the noise 
spectrum, and F is the filter function corresponding to 



VII. SUMMARY 

We investigated the protection of dissipative quan- 
tum state preparation processes against qubit dephasing 
by interlacing the DD control pulse sequences with the 
preparation control. The basic idea is to average out 
the noise channels over time, while preserving the irre- 
versible dynamics for generating desired entanglement in 
the steady state. By utilizing DD control consisting of 
sequences of short 7r-pulses, the dephasing noise can be 
suppressed to certain order of pulse interval and thus 
high-fidelity state preparation can be realized when the 
DD control has sufficiently small pulse interval. With the 
help of generalized Magnus series, we investigated the or- 
der of suppression of noise channels. For inhomogeneous 
dephasing, the leading noise term in the Magnus expan- 
sion is of second order of the pulse interval, and is largely 
determined by the average pulse interval rather than the 
types of DD control sequences used. The DD protec- 
tion efficiency is demonstrated by numerical simulations 
on two exemplary state preparation schemes for realizing 
many-body singlets and linear cluster states respectively. 
DD protection from dynamical noise is discussed for the 
example of Gaussian noise of a semiclassical description, 
where the leading noise effect is also of second order of 
the pulse interval. 
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